Evidence of a higher-order singularity in dense short-ranged attractive colloids 
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We study a model in which particles interact through a hard-core repulsion complemented by a 
short-ranged attractive potential, of the kind found in colloidal suspensions. Combining theoretical 
and numerical work we locate the line of higher-order glass transition singularities and its end-point 
— named A4 — on the fluid-glass line. Close to the A4, point, we detect logarithmic decay of density 
correlations and sub linear power-law increase of the mean square displacement, for time intervals 
up to four order of magnitudes. We establish the presence of the A4 singularity by studying how 
the range of the potential affects the time-window where anomalous dynamics is observed. 



Recently, a great deal of interest has grown around dy- 
namical phenomena arising in dense coUoidal suspensions 
when particles interact through a hard-core repulsion fol- 
lowed by a rather short-ranged attractive potential [1-3]. 
Some of these unusual phenomena, initially predicted 
theoretically [4-6], have been observed experimentally 
[7-10] and in numerical simulations [11-14,23]. The nov- 
elty, as compared to the well-studied hard-sphere-colloids 
case, comes from the possibility of generating structurally 
arrested states with an additional physical mechanism, 
driven by the range of the attractive interaction. Indeed, 
in the case of hard-sphere, structural arrest is controlled 
by the well-known excluded volume cage effect, where 
particle motion is hindered by the presence of neighbor- 
ing ones. This mechanism generates a localization length 
of the order of 10 % of the particle diameter a. In the 
additional mechanism, the motion of the particles is re- 
stricted by the adhesiveness, i.e. the formation of non- 
permanent bonds due to the attractive part of the inter- 
particle potential; particles are confined by bonds and 
the corresponding localization length is fixed by the at- 
traction range A. When A <C cr, the interplay of the two 
localization mechanisms creates liquid states for pack- 
ing fractions higher than those possible for a pure hard- 
sphere system. This means that the attraction stabilizes 
the liquid and leads to a reentry phenomenon, where this 
new liquid state can arrest into a glass by cooling as 
well as heating. The experimental and numerical confir- 
mation of this prediction [9-14] suggests that the mode- 
coupling theory (MCT) [15], which was the basis of the 
cited theoretical work [4-6] , can contribute to explain the 
slow dynamics in short-ranged attractive colloids. 

These systems are characterized by three control pa- 
rameters — the packing fraction cf), the ratio 9 of the 
thermal energy kgT to the typical well depth and the 
range A of the attractive potential. Within MCT, the 
phase diagram of this three-dimensional control parame- 
ter space (shown in Fig.l) is organized around a critical 
point {(p* ,6* , A*), referred to as a type A4 higher-order 



glass transition singularity in MCT classification. A4 is 
the end-point of a line of higher-order singularities (of 
type A3). From a physical point of view, A4 is char- 
acterized for being the only singularity point accessible 
from the liquid phase. For A > A* no singular points 
are predicted by the theory, while for A < A* the ^3 
singularity points are buried in the glass phase and their 
presence can be only observed indirectly [7,9,11,14]. Near 
the A4 singularity, MCT predicts a structural relaxation 
dynamics, which is utterly different from that known for 
conventional glass-forming liquids. It is ruled by loga- 
rithmic variations of correlators as function of time and 
subdiffusive increase of the mean squared displacement 
(MSD). Theory makes precise predictions for the time- 
interval where this unusual dynamics is expected, as well 
as for its variation with changes of the 4>, 9 and A pa- 
rameters. Indications of logarithmic dynamics have been 
reported [7,9,11,14], but neither a systematic study of 
the dependence on the control parameters has been pub- 
lished nor the A4 point has been located numerically. It 
is in fact difficult to discriminate between straight lines 
caused by some infiection points in the usual logarithmic 
representation of correlators and the logarithmic decay 
laws generated by the higher-order singularity. In this 
Letter, we establish the existence of the location of the 
A4 point for a well defined model, amenable of simultane- 
ous numerical and theoretical treatment. Reported data 
— with control parameters explicitly chosen close to the 
A4 point — exhibit the mentioned laws over time inter- 
vals up to four orders of magnitude. More importantly, 
we show that the decay patterns vary with changes of the 
control parameters and wave- vectors as expected [16,19], 
properly testing the theoretical predictions. 

We simulate a 50% — 50% binary mixture of = 700 
hard spheres of mass m, with diameters a^^ and (7^^ 
and ratio cr^^/a^g = 1.2. The hard core between parti- 
cles of different type a^^ = 0.5(0-^^ -I- a^^). The hard 
core potential is complemented by an attractive square 
well potential of depth Uo, independent on the particle 
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type [17]. The small asymmetry in the diameters is suf- 
ficient to prevent crystallization at high values of </>. We 
focus on two specific systems, which we label Si and S2, 
differing in the width of the square well Aij. Si and S2 
have Aij = O.OSlcTy, and Aij = 0.043c7ij respectively 
{i = A,B). The system has been extensively studied 
in previous simulations [14]. As discussed in the follow- 
ing, 5*2 is chosen to coincide with the critical amplitude 
parameter within the accuracy of our calculations. Av- 
erages over five independent realizations have been per- 
formed to reduce noise. Equilibration has been carefully 
checked. We note that our results refer to Newtonian 
dynamics, but arc relevant also to Brownian dynamics in 
the structural relaxation time window [21]. 




FIG. 1. Glass transition line for the system Si calculated 

through the vanishing of the diffusion (open symbols) and cal- 
culated with MCT(full line). The MCT ideal glass line for S2 
is given by the dashed line. The star indicates the location 
of the A4 singularity (on the 52 ideal glass line) while a filled 
square locates the point (on the Si line). The MCT curves 
have been transformed according to the mapping discussed in 
the text. The T and <j) of the simulated point arc indicated 
by the black dot. The inset shows the power-law fits for the 
isothermal diffusion coefficient D (data from Ref. [14]) for Si 
as a function of — </>c (see text). 

To estimate the predicted location of the A4 point for 
the binary mixture square well model (and consequently 
select the parameters to be used in the simulations) we 
proceed as follows. First, we calculate the glass line 4>c{T) 
for the mixture Si extrapolating diffusivity D data [14] 
according to D ~ (0 ^ 4'c{T))'^'^'^\ for eight different 
T. Second, we calculate the MCT ideal glass transi- 
tion line. The partial structure factors, the only input 
needed, are calculated with the Percus-Yevick approxi- 
mation, numerically solving the Ornstein-Zernike equa- 
tion [18]. Third, we map the theoretical glass line on 
the simulation line, following the procedure first used by 
Sperl [19]. Indeed MCT cannot reproduce the numerical 
values for (j)c. As shown in Fig.l, the linear transforma- 
tion (j) 1.897 (f> - 0.3922 and T 0.5882 T - 0.225, al- 



lows to superimpose the MCT result with the simulation 
data, in the studied region of T and (j). By calculating 
the MCT ideal glass transition line for several values of 
A, the theoretical location of the A4 point for the binary 
mixture square-well model is obtained. Assuming that, 
to a first approximation, the same linear transformation 
holds also for the critical value of the well width, we find 
that the A4 location maps to A*^- = 0.043crjj, (p* = 0.611 
and 9* = 0.416. This allows us to select a state point 
with parameters close to the critical ones and perform a 
simulation (in thermodynamic equilibrium) close to A4. 
We choose [20] 6 = 0.4 and (f) = 0.6075 and compare, for 
the same T and cj), the dynamics for A*^ (system S2) and 
for the close-by value A^j = 0.031(Tjj (Si). 




t 

FIG. 2. Mean squared displacement of the A particles, 
{5r\), showing sub-diffusive behaviour within the marked in- 
tervals. At long times, the diffusive behavior {5r\) ~ i is 
recovered. 

Fig. 2 shows the MSD for the larger species (fr^) for 
both systems. We notice the presence of a subdiffusive 
regime at intermediate times, i.e. a variation according 
to ln{{6r\)) = f + a ln{t/T) for about three decades for 
^i, where ai = 0.44. The logarithmic regime extends by 
more than a decade if the attraction range is increased to 
that for system ^2. Simultaneously the exponent a de- 
creases to a2 = 0.28. Fig. 3 shows the density-fluctuation 
auto-correlation functions ^q{t) for a representative set 
of wave vectors q. The decay curves in Fig. 3 do not 
show the two steps scenario with a plateau characteristic 
of conventional glass forming liquids. Indeed it is im- 
possible to fit these curves for large q with the standard 
stretched exponential function. Instead, there is a region 
of clear logarithmic decay at ql^gg = 23.5 for Si and at 
~ 16.8 for ^2. The time-intervals over which log- 
arithmic decay is observed are of similar size as those for 
the MSD. For q < q*, the $g vs In t curves are concave, 
and for q > q* convex. 
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To show that the above described features are consis- 
tent with the ones predicted by MCT, we cite the general 
asymptotic decay law for the correlation function of a 
generic variable X, near the higher-order singularity [16] 

^x{t) = fx-hx [B^^Hnit/r] + Bfln''{t/T)\ . (1) 

Here r abbreviates a time-scale which diverges if the state 
approaches the singularity. The formula is obtained by 
asymptotic solution of the MCT equations, using the pa- 
rameter differences 0* — 0, 6** — 0, A* — A as small quan- 
tities, say, of order e. The coefficient S^^^ is of order 
y/e while i?^^-* is of order e. The amplitude hx is in- 
dependent of e. The first term fx is the sum of the 
non-ergodicity parameter of variable X at the singular- 
ity and a correction of order e. Terms of order e^/^ are 
neglected [16]. We can interpret ln{5r\) also according 
to Eq. (1). Hence, the straight line for with slope a\ 
shown in Fig. 2 is consistent with the assumption that 
is close to the singularity. Changing from Si to 52, 
the slope a = hj^^j^B^^^ has to decrease and the range of 
validity of the leading order description has to expand. 
This is demonstrated impressively by the data. 
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FIG. 3. A-type particle density correlators for (top 
panel) and S2 (bottom panel), showing six different wave vec- 
tors (from top to bottom qa^g = 6.7, 11.7, 16.8, 23.5, 33.5 and 
50.3). For correlators indicated by an arrow, a logarithmic 
behavior is observed within the selected time window. The 
dashed lines show the fits according to Eq.l. The vertical 
dashed lines indicate the fitting interval. 



To estimate the possibility of describing study the time 
dependence of ^q{t) according to Eq. 1, we fit the den- 
sity autocorrelation functions to a quadratic polynomial 
in log{t/T) for different q values. Fits are reported in 
Fig. 3. The fitting time window extends from about 2.5 
decades for the state point Si to four decades for ^2. 
The fitting parameter shown in Fig. 4, provides an es- 
timate of the non-ergodicity parameter at the A4 point. 
We find that fq does not depend on the state point. This 



confirms the preceding conclusion that the studied state 
points are very close to the singularity and that the order 
e correction in fq cannot be detected. 
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FIG. 4. The fitting parameters fq, H^'^^ and H^^' of the 
asymptotic logarithmic law, for Si and 52 . The central panel 
also shows that a multiplication by 2.43 of Hq^-' for the 52 sys- 
tem gives Hq^^ for Si, confirming the factorization of Hq^\ 
At the wave vector where Hq^^ = (qJ^'bs = 23.5 for Si 
and qiugg = 16.8 for ^2) correlation functions display a pure 
logarithmic decay. 
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The fit parameters for the coefficient Hq 
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hqB^^) 

are reported in the middle panel of Fig. 4. The decrease 

of h'^^ upon changing from 5*1 to 5*2 is in agreement 
with the prediction that S^^' tends to zero with e ap- 
proaching zero. As shown in Eq.l, factorizes into a 
control-parameter independent factor hq, depending on 
q, and a control-parameter dependent factor B^^', inde- 
pendent of q. This implies that the g-dependence of Hq 
should be the same for and ^2. As shown in Fig. 4, this 
property is verified by the data. The same property docs 
not apply to Hq ' = hqBq ' because of the g-dependence 



(1) 



of B^K Moreover, the fitting results confirm that Hq 
is smaller than Hq'\ as expected being the first of or- 
der e and the second of order ^/e. The wave- vector value 
(2) 

where Bg ' = allows to identify the characteristic length 
scale associated to the pure logarithmic decay (see Fig. 3). 
Such length scales are much shorter than the typical first 
neighbor shell. 

To explicitly test the mentioned factorization, Fig. 5 
shows a set of correlators rescaled to ^q{t) = {^q{t) — 
fq)/H^q^\ MCT predicts that in leading order e all 4,(t) 
vs ln{t) curves collapse on the logarithmic decay law 
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—ln{t/T). This is indeed shown for a large time inter- 
val that expands upon approaching the singularity, as it 
is demonstrated by comparing results for S\ and ^2 [22]. 



minating discussions, G. Foffi and M. Sperl for comments 
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FIG. 5. Scaled correlation functions $q(t) (see text) vs 
ln{t/T) for the two samples and selected values of qas around 

(2) 

the characteristic values q* for which Hg is close to zero. The 
vertical dashed lines indicate the time interval over which the 
scaling is observed. The values of t are 10 and 10^ for Si and 
S2 respectively. 

Results reported in this Letter have an intrinsic value 
associated to the observation, in a particularly simple 
system — a square well potential — , of a particularly 
complex dynamics over more than four decades in time. 
They show that MCT, a theory essentially developed to 
address the problem of the excluded volume glass tran- 
sition is able — without any modification — to handle 
the logarithmic dynamics and to provide an interpreta- 
tive scheme in term of A4 point . The numerical results 
and the comparison with the theoretical predictions do 
constitute, in fact, a stringent verification that the loga- 
rithmic dynamics in the density autocorrelation functions 
and the subdifi^usive behavior in the MSD can be fully ra- 
tionalized by MCT. Indeed, data in Figs. 3 and 4 are in 
agreement with MCT predictions (Fig. 12 and 13 of Ref. 
[19]). The simplicity of the model and the complexity of 
the dynamics suggest that the short-range attractive col- 
loids have the potentiality to become a benchmark for the 
development of extended theories of the glass transition. 
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